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Abstract 

The pricing of options, warrants and other derivative securities is one of the 
q ■ great success of financial economics. These financial products can be modeled 

and simulated using quantum mechanical instruments based on a Hamiltonian 
£Sj \ formulation. We show here some applications of these methods for various 

potentials, which we have simulated via lattice Langevin and Monte Carlo 
algorithms, to the pricing of options. We focus on barrier or path dependent 
options, showing in some detail the computational strategies involved. 
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1 Introduction 



Since Black and Scholes' (BS) option pricing model gained an almost immediate ac- 
ceptance among the professional and academic communities, the trading of derivative 
securities skyrocketed. Derivative securities (such as options) are financial securities 
whose payoffs depends on other underlying securities, and the BS model was the first 
universally accepted modeling of these financial instruments. 

However, in recent years, financial engineeers 1 have created a variety of complex 
options that are collectively called exotic options. 

The payoffs on these options are considerably more diverse than the payoffs on 
standard BS options or on other straightforward generalizations of them. Some of 
these financial instruments, widely used in complex portfolios incorporating thou- 
sands of these (correlated) elementary instruments, may be easier to analize using 
standard quantum mechanics. 

While we focus on the possibility of applying fundamental theories to the real 
economy, we have to remark that the pattern seems to be much wider than expected, 
since even quantum field theory methods [T] have found their way to the financial 
world. 

Most of the mathematical methods involved in the analysis of financial systems 
have been based, so far, on the simulation of stochastic processes by diffusion equa- 
tions coupled to stochastic sources, i.e. stochastic equations of Langevin type. More 
recently, there has been an interest in the analysis of various financial instruments 
using the path integral formulation. 

Use of a path integral formulation has some advantages. First, it is in close rela- 
tion to the lagrangean description of diffusion processes, second, it opens the way to 
the use of quantum mechanical methods, on which we briefly elaborate. After a de- 
scription of the path integral in the Black Scholes model, which we study numerically, 
we turn our attention to the analysis of barrier options. Barrier options are studied 
here by simulating an artificial quantum mechanical model in which a potential V(x) 
is added to the Black Scholes lagrangean, as first suggested in ref. fT]. Simulations 
are carried out using both Langevin and Monte Carlo methods, comparing in some 
limiting cases, where possible, our numerical results with analytical ones. 



2 Langevin Evolution 

In the top down description of theoretical finance, a security S(t) follows a random 
walk described by a Ito-Weiner process (or Langevin equation) as 

= 4>dt + aR{t)dt, (1) 



1 Under this professional name perform actual research a large number of physicists and mathe- 
maticians, with a remarkable number of high energy physicists 
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where R(t) is a Gaussian white noise with zero mean and uncorrelated values at time 
t and t' (R(t)R(t')) = 5(t — t'). (f> is the drift term or expected return, while a is a 
constant factor multiplying the random source R(t), termed volatility. 

As a consequence of Ito calculus, differentials of functions of random variables, 
say f(S,t), do not satisfy Leibnitz's rule, and for a Ito-Weiner process with drift (0) 
one easily obtains for the time derivative of f{S,t) 

d f df 1 2Q2 d 2 f df df 

Tt = m + 2 aS dsi + (l)S ds + aS ds R - (2) 

The Black-Scholes model is obtained by removing the randomness of the stochastic 
process shown above by introducing a random process correlated to (j2J)- This oper- 
ation, termed hedging, allows to remove the dependence on the white noise function 
R(t), by constructing a portfolio II, whose evolution is given by the short-term risk 
free interest rate r 

~dt = rIL < 3 > 

A possibility is to choose IT = / — ^S. This is a portfolio in which the investor holds 
an option / and short sells 2 an amount of the underlying security S proportional to 
|£. A combination of (j2J and (jHJ) yields the Black-Scholes equation 

There are some assumptions underlying this result. We have assumed absence of 
arbitrage, constant spot rate r, continuous balance of the portfolio, no transaction 
costs and infinite divisibility of the stock. 

The quantum mechanical version of this equation is obtained by a change of 
variable S = e x , with x a real variable. This yields 

with an Hamiltonian Hbs given by 

Notice that one can introduce a quantum mechanical formalism and interpret the 
option price as a ket |/) in the basis of \x), the underlying security price. Using 
Dirac notation, we can formally reinterpret f(x,t) = (x\f(t)), as a projection of an 
abstract quantum state \f(t)) on the chosen basis. 

In this notation, the evolution of the option price can be formally written as 
\f,t) = e tH \f,0), for an appropriate Hamiltonian H. 

2 short selling of the stock should be possible 
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3 Options and Barrier Options 



3.1 Generalities 

Let the price at time t of a security be S(t). A specific good can be traded at time t 
at the price S(t) between a buyer and a seller. The seller (short position) agrees to 
sell the goods to the buyer (long position) at some time T in the future at a price 
F(t,T) (the contract price). Notice that contract prices have a 2-time dependence 
(actual time t and maturity time T). Their difference r = T — t is usually called 
time to maturity. Equivalently, the actual price of the contract is determined by the 
prevailing actual prices and interest rates and by the time to maturity. 

Entering into a forward contract requires no money, and the value of the contract 
for long position holders and strong position holders at maturity T will be 

(-iy(S(T)-F(t,T)) (7) 

where p = for long positions and p — 1 for short positions. Futures Contracts are 
similar, except that the after the contract is entered, any changes in the market value 
of the contract are settled by the parties. Hence, the cashflows occur all the way to 
expiry unlike in the case of the forward where only one cashflow occurs. They are 
also highly regulated and involve a third party (a clearing house). Forward, futures 
contracts and, as we will see, options go under the name of derivative products, since 
their contract price F(t, T) depend on the value of the underlying security S(T). 

In the simplest option, such as a call option, we have seen that the payoff function 
is defined to be the value of the option at maturity time (r = 0). Therefore, the 
specific path followed by the underlying security is not relevant in order to establish 
the price at maturity, except for its final value. 

Barrier options are, instead, path-dependent. This means that the payoff is de- 
pendent on the realized asset path, and certain aspects of the contract are triggered 
if the asset price, from start to end of the contract, becomes too high or too low. 

Barrier options are very popular for various reasons. An investor may have very 
precise views about the behaviour of a security or he may use them for hedging 
specific cashflows, to decide to purchase them. In the following, when comparing 
path dependent options to the simplest options, such as standard calls or puts, we 
will refer to the latter as to vanilla options, using a common financial jargon. 

3.2 Terminology and Definitions 

There are some advantages -and natural limitations- in purchasing a financial instru- 
ment such as a barrier option. If the purchaser wants the same payoff typical of a 
vanilla option, but believes that the upward movement of the underlying will not be 
likely, then he may decide to buy an up-and-out call option. The cost of this contract 
will be cheaper than the purchase of a corresponding plain vanilla option, but there 
will be severe limitations on the upward movement of the option. 
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The physical picture of an up-and-out option is that of a brownian motion of the 
underlying asset (x) that is immediately killed as soon as the asset hits (from below) 
the barrier B (x = B), which is specified in the contract. 

Similarly, a down and out provision renders the option worthless as soon as the 
asset price hits a barrier B from above. The payoffs in the two cases are given by 

guo( x ,K) = max(Sr — K)8(B — x) 

g DO (x,K) = max(S T - K)8(x - B) (8) 

for a up-and-out (UO) and a down-and-out (DO) option call respectively. Here, 9() 
denotes the standard step function. A terminology used to describe contracts with 
these features is knocked out options. In contracts of this type it is agreed there will 
not be any payoff if the barrier B is hit. 

Similarly, the market offers contract with additional limitations on the allowed 
variation of the underlying asset. For instance, double knock out options have restric- 
tions on the asset variability delimited by two barriers (B\ < B u ) both from above 
{By) and from below (B{), and give zero payoff if any of the two barriers is hit by 
the asset from inception time t to expiry time T. 

Knock in options are dual, in an obvious sense, to knock out options. Knock in 
options, in fact, are contracts that pay off as long as the barrier B is hit before expiry. 
If the barrier is hit, then the option is said to have knocked in, otherwise their payoff 
is null. 

Furtherly categorizing these latter types of options, the position of the barrier 
respect to the initial value of the underlying allows to distinguish between up-and-in 
options and down-and-in options. The payoffs of these contracts are given by 

gui(x, K) = max^Sx — K)9{x — B) 

g DI (x,K) = max(S T - K)6(B - x). (9) 

For definiteness, in the analysis that follows up, we will focus our attention to 
knocked out payoffs of the types described in eq. (jEJ). 

In knocked out options, single or double, killing of the brownian motion is, needless 
to say, instantaneous, and takes place as soon as the brownian motion of the asset 
hits any of the barriers. 

This aspect of the contract is an unpleasent feature since it introduces a disconti- 
nuity in the dynamics, with attached risk management problems both for option buy- 
ers and sellers. Such risks, for instance, are those due to erroneous price movements, 
or to an instantaneous spiky behaviour of an asset, moving upward or downward and 
penetrating a given barrier, which can lead an investor to the loss of all his invest- 
ment. In other unpleasent situations, when large positions of options accumulate in 
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the market and are all characterized by the same barrier, trading can drive the asset 
to the barrier, generating massive losses. 

There are various ways by which more conservative and safer contracts can be 
defined, while maintaining some of the features of knock out options. This is achieved 
by introducing a finite knock out rate, thereby smoothing out the effect of the barrier. 
Our goal is to show how it can be implemented in a self-consistent path integral 
formulation and characterize the pricing of these path dependent options. 



4 Quantum Methods in Finance 

To establish a path integral description of a stochastic process we need a lagrangean 
and the corresponding action. This can be easily worked out for the BS model, 
starting from the Hamiltonian given in eq. (|B|) . We easily gets 



1 / dx 1 



'BS 



2 



2^{ii +r -r)- r < 10) 

and the corresponding action, expressed in terms of time to maturity r 

S BS = [ L BS (f)dt' (11) 
Jo 

which can be used to define a corresponding path integral for a fictitious quantum 
mechanical process in the variable x, the logarithm of the underlying asset 



/+oo 
d: 
-oo 



x(t)e s[x] (12) 

with the boundary conditions x(ti) = Xi and x(tf) = Xf. The variable x = log(S') 
which identifies the quantum mechanical state of the system will be refered to as to 
the stock price. The pricing kernel for the stock price is given by the 



p BS (x,x',r) = J DX BS e SBS 
= (x\e- rHBS \x') 



(13) 

with 

r rex 
DX BS = U T t=0 / dx(t). (14) 
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4.1 Generalized Potential 

For barrier options it is tempting pp to introduce a potential V(x) in order to set up 
a constraint on the stochastic process described by the stock price x. 
The corresponding generalized Hamiltonian now reads 

H v = + (\a 2 - V(x)) |- + V(x). (15) 



2 dx 2 \2 / dx 

This Hamiltonian is equivalent to a stochastic process, as given in ([86)1 . with dis- 
counting done by exp(— J dtV(x(t))). 

It can be shown pQ that Hy obeys the martingale condition, and hence can be 
used for studying processes in finance. 

The non-Hermiticity of Hy is of a particularly simple nature, and it can be shown 
PP that for arbitrary V, Hy is equivalent by a similarity transformation to a Hermetian 
Hamiltonian Hes 3 given by 

(16) 

+ i^ + J-y 2 +Iy+^ (it) 

2 dx 2<t 2 2 8 1 ; 

rfyF(y) (18) 

Note that i^ES is Hermetian and hence its eigenfunctions form a complete basis; 
from this it follows that the Hamiltonian Hy can also be diagonalized using the 
eigenfunctions of i^ES- In particular 

H ES \4> n > = E n \4> n > (19) 

^H v \^ n > = E n \^ n > (20) 
where 

|^n> = e s \<p n > (21) 

<ip n \ = e' s < (f) n \ ^< ip n \ (22) 

For the Black-Scholes Hamiltonian Hbs we have V(x) = r and hence 

H BS = e s H ES e~ s (23) 

2 dx 



Hes = 


e 




where 






Hes = 




a 2 <9 2 




2 &r 2 


and 








1 


1 


s = 


2 


x ~ — 
a 1 



(XX 



e-[- V ^+7]e- M (24) 



where 



7 = 2a 2 ; a = ^*— (25) 



3 Note that for more complex Hamiltonians such as the Merton-Garman Hmg finding the equiv- 
alent Hermetian Hamiltonian Hes is far from obvious 
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5 Path Dependent Options 



We can have very complicated path dependent options since an option is an arbitrary 
random variable on the underlying sample space, or in other words, a completely 
arbitrary functional of the history of asset prices. For many but not all kinds of path 
dependent options, we can extend the technique of obtaining from the path integral 
a Hamiltonian for a quantity related to the option which is path independent and 
can therefore be represented as a wave function. The solution for the pricing kernel 
of this quantity then gives us the solution for the path dependent option. It must be 
noted that this quantity cannot be a traded asset as all traded assets evolve with the 
Black-Sholes Hamiltonian 

Hbs = ~ rS ds - —dsi (26) 

according to 

-*+!)/ = " <») 

in the Black-Scholes model. Let us first look at how this is done for some relatively 
simple (but more complicated than simple barrier options) path dependent options. 

5.1 Soft barrier options 

These options have been considered in detail in Linetsky jS] . They are similar to the 
barrier options considered above but do not knock out the option completely when 
the barrier is hit. Instead, for soft barriers one discounts the final payoff, at some 
rate, by the exponential of the amount of time spent inside or outside the barrier. 
For example, for a down and out barrier is a discounted step option whose barrier is 
at B and strike price at K, the payoff at expiry is 

- Vtb -(S t -K) + (28) 



e 



where tb_ is the time spent below the barrier B and V is the discounting factor. 
Considered as a path integral, the current price of the option is given by 



Defining a potential 



r-x(T)=x' 

dx'(e x ' - K) + / Vxe SBS e~ VTB - (29) 

Jx(0)=lnS 



• V x < In B , , 

V{X) = <0 *>>n B (30) 



we see that that the path integral is equivalent to 

•x(T)=x' 



dx'(e x ' - K)+ I Vxe s (31) 

x(0)=lnS 
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with the action S now being given by 



S = 



JdtL{x,x) = -^Jdt^(x-r-^j + r + V(x)j (32) 



In other words, we have just introduced a potential into the problem. The Lagrangian 
is now 

L = L BS + V(x) (33) 
and the Hamiltonian now has an extra term — V(x). The new Hamiltonian is therefore 

a 2 d 2 fa 2 \ d 

^-¥^+(~'j&+ r + V ' W (34) 

Hence, the solution of the step option price in the Black-Scholes model is equivalent 
to the solution of the plain vanilla option in a model with the above Hamiltonian. 
More generally, the pricing of an option whose final payoff is 

e-S v ^ dt (S T -K) + (35) 

is equivalent to the pricing of a plain vanilla call option in a model where the Hamil- 
tonian is 

H = H BS + V(x) (36) 

If we can find the eigenvalues and eigenfunctions of the operator H, we can write 
down the pricing kernel using the decomposition 

(x\e- T "\x') = ^2(x\n)(n\e- T(l \n)(n\x') = ^ e^"^')^) (37) 



where the eignevalues of H are E n and the eigenfunctions corresponding to these 
eigenvalues are ip n . When the eignenvalues are continuous, the sum becomes an 
integral as is the case in the calculations for the single barrier options. In this case, 
we can consider the Laplace transform of the pricing kernel 

'ST / I -tH\I 



dTe- ST (x\e- rH \x'} (38) 



o 



which is seen to be the Green's function of the operator s + H (this can also be 
directly seen from the Feynman representation of the pricing kernel). Once we find 
the Green's function, we can perform the inverse Laplace transform to get the pricing 
kernel and hence the solution to the problem. 

To see how this works, let us take the example of the simple barrier option. We 
modify the state space using the terms a and (5 so that we only deal with a standard 



8 



Brownian motion. In that case, the potential only influences the boundary conditions, 
so we have to find the solution to the equation 

1 d 2 



o _ -s)G(x,x';s) = -5(x-x') (39) 

with the boundary conditions G(x,x';s) = 0,x,x' = b = \nB and 
linXr-Kxj G(x, x'; s) = 0. The Green's functions can be easily found using standard 
methods and the result is given by 



„, , x 2sinh v / 2i(a;-&)e- v ^( 1 '- 6 ) , . 
G{x,x ; s) = -= B(x — x) 

v 2s ^ (40) 
2 sinh y/2a(x! - 6) e -v^^-«>) 
H v _ ' Q(x - x') 

whose inverse Laplace transform is the pricing kernel 

p BS (x, r; x') - -_L= e -^-a(*-*') exp [ - -^(x + x' - 2B) 2 ] (41) 

where pbs is the Black-Scholes pricing kernel and where the adjustment for the 
transformation has been made. 

This technique is applied to find a closed form solution for the step option in 
Linetsky [Oj. If we choose the variables such that the stock price S = Be ax so that 
the barrier is at x = and again only deal with the standard Brownian motion, we 
find that the Green's function is given by 

1 ( e VTs\x-x'\ _ V 's+V -y^ c ^Ts(x+x')\ X X' > 

V2S \ ^S + V + y^S J 



G(x, x'\ s) 



3 y/2(s + V)x-V^x' 
y /2(s+V ) + y/2^ 



x < 0,x' > 
x > 0,x' < 



^2{s+V)+VTs 

' y/2{s+V)\x-x'\ _ V^W~^ c ^/2(s+V)(x+x')\ % % , < q 



^ y/2(a+V) 

whose inverse Laplace transform gives us the pricing kernel. 



(42) 



5.2 Asian options 

There are several options which cannot be put into a simple form by discounting 
alone. One such option which is also fairly popular in the market is the Asian option 
which has been considered in detail in the literature. The Laplace transform of an 
out of the money Asian option was found in Geman and Yor [4 . The payoff of the 
Asian option is defined to be 

max(o,^J S(t)dt-K) (43) 
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We can write the option price as a path integral 

J dp J dv J Vxe SBS e tp{ ± f« dteX(t) - u \v - K) + (44) 

using a standard expression for the Dirac delta function. We see that we can consider 
this expression as a plain vanilla call option with a Lagrangian modified by — % ^e x . If 
it were not for the term i, this would be a reducible to a relatively standard problem. 
However, since the additional potential term is now complex, the solution is not so 
simple. 

One option which is a fairly good approximation for the Asian option but which 
is easily solvable is the geometric Asian option. Its final payoff is defined to be 

max(o,e^o - K^j (45) 
To solve this, let us write the Black-Scholes evolution as 

% = (r-C)+*v(t) (46) 



dt V 2 
where r)(t) is white noise. Hence, 



x(t) = x(0) + (r - y ^ t + jf dt'r](t') (47) 



(48) 



and 

- J dtx(t) = x(0) + - [r - °- J + £ jf (T - i)77(*)* 

To find the distribution of the last term, we make use of the generating function for 
white noise to give 

— / Vrje'^ £ dtr, 2 (t) e ip(% /„ T ' dt(T-t) V (t)-u) ^ 

2-k J 1 

which evaluates to 

I 3 3. 2 , 

V^r e ~^ (50) 

In other words, the distribution of ^ x(t)dt is N(0, ^y-). Therefore, the price of 
the geometric Asian option is given by 

c = SW(di) - Ke- r{T - t] N(d 2 ) (51) 

where 

V3(^(#)+(i + f|)(r-0) V3 (in (I) + (l - £ ) (T - t)) 

di = — ^ ; ' '-, d 2 = — ^ ) J J - 52 
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5.3 Seasoned options 

In the above discussion, we have only discussed how to value path dependent op- 
tions where the path dependence starts at the present. In practice, the problem of 
solving for the price of path dependent options after they have been initiated is very 
important. Such options are called seasoned options. In many cases, the valuation 
of seasoned options proceeds very similarly to that of new options. 

Let us consider a seasoned soft barrier option with discounting by a potential 
V(x) at time t > where the path dependence has started at zero time. We denote 
the maturity time be T. The value of this option is given by 

.-S^j^C*™^-^ (53) 

where we should take into account the discounting until time t separately since the 
history is already known. We see that apart from this factor, there is no substantial 
difference between the valuation of the new and seasoned options. 

For Asian options, we see that we can value the seasoned option provided we can 
value the new option since the probability distribution of the remaining part of the 
average will determine the option price at the time when combined with information 
about the contribution to the average of the revealed historical price. Hence, we see 
that we can price seasoned Asian options if we can price new Asian options. 



6 Solving the double knock out barrier option 

A double barrier option is an option whose value reduces to zero whenever the price 
of the underlying instrument hits the barriers which we denote by e a and e b . Hence, 
the price of a double knock out barrier European call option expiring at time T and 
with strike price K at time to provided it has not already been knocked out will be 
given by 

e- r ^E t [(e^ - K) + ]l a<x{t)<bM<tl<T (54) 

where 1 stands for the indicator function. It is sufficient to solve for the probability 
distribution of x(T) for those paths which do not go outside the barriers (in other 
words, the pricing kernel). 

Written as a path integral, the formula is 

.-^./^-.W-^^-JO, (55, 
where Sbs is the Black-Scholes action 

S bs = ~^- 2 J dt(x + r-^f (56) 
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While the step functions look complicated in the path integral, they can be seen 
to be having the effect of an infinite potential barrier since they effectively prohibit 
the path from entering the forbidden region outside the barriers. Hence, the problem 
might be better solved using the Hamiltonian and this is indeed the case. 

In the Schrodinger formulation, the above problem is to find the pricing kernel 
for a system with the Hamiltonian 

H = H BS + V{x) (57) 

where the Black-Scholes Hamiltonian is given by 

°" 2 9 .a 2 . d 

*- = -Ta^ + <Y-'% < 58 > 

and the potential V(x) is given by 

V(x) = { a < x < b (59) 




This is very similar to the well known problem of a particle in an infinite potential 
well except that the Hamiltonian has an extra term involving Jj^ which makes it non- 
Hermitian. 

This problem can be solved by transforming the underlying wave functions. By 
making the transformation (x | 0) = e~ a ( x ~ aS) (x | t/j) and ((f) | x) = e a< - x ~ a ^(ip | x), 
where |0) are the vectors in the new (Hilbert) space, and (^l are the original 
vectors and their duals respectively and a = a ^1~ r . In this new space, the Black- 
Scholes Hamiltonian takes the simple Hermitian form 

The problem is now identical to that of a quantum mechanical particle of mass 
-js (in units where h — 1) in an infinite potential well. As is well known in this case, 
the allowed momenta are p n = The eigenfunctions are hence given by 



(x\^ n ) = e ai - x ~ a \x | 4> n ) = \\^- a ^ a( - x ~ a) sin p n (x - a) (60) 



$ n \x) =e-^ x ~ a \<p n \x) =-J-^—ie- a ^- a hmp n (x-a) (61) 

V o — a 

where (x \ <p n ) are the eigenfunctions of the quantum mechanical particle in an infinite 
potential well. 
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The eigenfunctions are orthonormal and form a complete basis since 
( x I | x') = -—e a(x '" x ' ) ^2smp n (x - a) sinp n (x' - a) 

1 O 

n=l n=l 

1 cy(r-r') ( ^ n1X I l\ ^ n7T / I n \ 

= 5(F^ e< £ (exp— (z-zj-exp— (z + z -2a) ; 

v ' n=— oo v 

= J0_ e a(x-x>) ( 5 ( n(x-x') \ _J 7r(x + x' -2a) \ ' 
b— a \ \ b—a J \ b—a J 

= S(x - x') 

since a < x < b and a < x' < b. 

The pricing kernel is hence given by 



oo oo 



(x\e- TH \x')^J2( x \^)^n\e- TH \M(^n'\x') 



n=l n'=l 

oo 



Y J {x\i>n){i>n\x')e- TE " 



n=l 



1 / ra 2 

exp I h a(x — x t 



2(6 -a) 

^ exp r °" 2 P» VgipnCx-xQ _ e ip„(x+s'-2a)) 
n=— oo ^ ' 



iyisix — x ) iynix + x — 2a) 

exp exp 

o — a ' b — a 



where 



v - / (x - + 2n(b - a)) 2 (x + x'-2a-2n(b-a)) 2 

z2 H w exp 2^ 

(63) 



, = (64) 



and the identity 

5{y -n)= ^ ( 65 ) 



oo 
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has been used. 

Hence, we see that the pricing kernel (apart from the drift terms) is given by an 
infinite sum of Gaussians. To check its reasonableness, we check the value in the 
limits b — > oo and a — > — oo. In the former case, only the n = term contributs and 
in the latter, only the n = and n = 1 terms contribute. It is easy to see that, in 
both cases, the result reduces to the solution for the single knockout barrier pricing 
kernel. When both limits are simultaneously active, only the first term in the n = 
term exists and it is easily seen that gives rise to the well known Black-Scholes pricing 
kernel. 

We can now evaluate the price of a double barrier European call option using the 
pricing kernel from The result is seen to be 



f = J2 {e~ 2na{b ~ a) (e 2n ^ a) SN(d nl ) - Ke- rT N(d n2 )) 

n=—oo 

/ 2a 

- s 2a e - 2a{n{b - a) ~ a) ie 2n{b - a) —N(d n3 ) - Ke- rr N(d n4 ) 



(66) 



where 



ln(f ) + 2n(b-a) + r (r + ^) 
d nl = = ^ J - (67) 



(JWT 



l-n2 



^n3 



ln(f ) + 2n(b-a) + r (r - ^) 



CT\/T 



ln(g) + 2n{b - a) + r (r + 
ln(f|) + 2n(b _ a) + r (r - f 



CTa/T 



dni - (Ty/r (68) 

(69) 

= d n3 - (70) 



7 Monte Carlo Simulations 

The pricing kernel is the fundamental quantity to compute using the functional inte- 
gral. Related attempts can be found in the literature [3] . We assume a discretization 
of the time to maturity r in intervals e = t/N, with N an arbitrary (large) integer. 
For instance, for the Black Scholes model one gets the action 

N 

S BS = eJ2LBs(*) (71) 

i=l 

with 

WlH __l(£i^zI + r _^ 2 (72) 
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where we have introduced discretized positions for the variable 

For this purposes, we have used a standard Metropolis algorithm. If thermaliza- 

tion is slow, it is possible to resort to use sequentially Metropolis updates and cluster 

updates. The latter is an update for the embedded Ising dynamics in the lattice 

variables Xi/\xi\ (Swendsen-Wang, Wolff), and is included in for a faster generation 

of the thermalized paths of the stock price x(t). 

For processes involving a stochastic volatility (y = \og(V)) the expression of the 

path integral is more complicated and can be found in [2j. From now on we will just 

consider the case of a constant volatility. 

If we denote by g(x, K) the payoff function, with a strike price K, in this case the 

value of the option (its price) is given by the Feynman-Kac formula 



In actual simulations, it is convenient to compute directly the option price rather 
than the propagator itself. The simulation is done by taking the initial point x fixed, 
and letting the final point evolve according to its quantum dynamics. In this way 
a path (x, x') is generated. After the first thermalization, x' is allowed to undergo 
quantum fluctuations, at fixed x. Each x' is then convoluted with the payoff function 
and an average is performed. Finally, this procedure is repeated for several x values, 
so to obtaint the option price at time to maturity r. 

Figs. 1, 2 and 3, illustrate some simple results obtained by the monte carlo 
method. For illustrative purposes, we show the behaviour of the Black-Scholes model. 
Fig. 1 shows a typical thermalized path, generated from a given initial value x (at 
current time t = r) assuming a maturity of 300 days, while in Fig. 2 we have plotted 
several path for different starting values x of the stock at current time r. We have 
chosen an interest rate r = 0.05 and a 12 percent volatility a. Finally, in Fig. 3 we 
compare the analytical and the numerical evaluation of the Black-Scholes option price 
with a low resolution for ([73)1 . in order to separate the two curves, which otherwise 
would overlap completely, in order to illustrate the convergence of the Metropolis 
algorithm. 

8 Langevin Calculation of Option Prices with Po- 
tential: Numerical Methods 

As is well known from quantum mechanics, the path integral 



f(t, x) = dx'( 



^ HBS \x')g{x',K). 



(73) 



J — oo 




(74) 



with the Lagrangian L given by 




V(x) 



(75) 
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Figure 1: An example of thermalized path obtained from the simulation of the path 
integral (Black-Scholes) with r=0.05 and a = 0.12 



is the pricing kernel for a one-dimensional system of a particle of mass m moving 
in a potential V(x). This is very similar to the path integral for standard Brownian 
motion W{t) given by 

r VWe-y W2 ^ dt (76) 



The path integral for a more general stochastic process described by the stochastic 
differential equation (in physics, the equation is usually written out with everything 
divided by dt and with dW(t)/dt replaced by rj which represents white noise and the 
equation is then called the Langevin equation) 



dx(t) = a(x)dt + a(x)dW(t) 



(77) 



The Lagrangian can be found by solving the above for W . One easily obtains the 
path integral 

" Vxe -J dt l^(*-^) 2 (78) 

When we perform a Wick rotation for the action in (j74j) . we get the Lagrangian for 
Brownian motion if V(x) = 0. Hence, a free quantum mechanical particle can in 
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Figure 2: Several thermalized paths for (Black-Scholes) with r=0.05 and a = 0.12 

some sense be considered to be undergoing Brownian motion and can be modelled 
by the stochastic differential equation 

dx(t) = \ —dW(t) (79) 
V m 

The path integral for the Euclidean action can be numerically simulated using the 
Metropolis algorithm. Another method is to use the analogy and directly integrate 
the stochastic differential equation (|79|). The latter is more efficient as one does 
not have the problem of correlation between succesive configurations which reduces 
the accuracy of Monte Carlo calculations. There does seem to be a problem in 
that potentials cannot be included. However, this can be handled by including a 
killing term in the stochastic differential equation whose connection with the potential 
becomes clear when we consider the Hamiltonian. The integration of stochastic 
differential equations can also give efficient numerical calculations in quantum field 
theory, especially for gauge invariant theories where gauge fixing is trivial in this 
framework. This is dealt with in great detail in Namiki 

In the case of stock option pricing, we are usually dealing with the Black-Scholes 
stochastic differential equation 

dS(t) = rS(t)dt + aS(t)dW(t) (80) 
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Figure 3: Call option price for strike price 3 versus the logarithm of the initial value of 
the stock xo = log(So). the parameters are fixed as in figs 1. Shown is the analytical 
result vs the monte carlo result, with a low resolution of 10,000 configurations 



whose Hamiltonian is given by 



and Lagrangian by 



1 2 Q 2 d 2 Q d 
2 ab dS^ + rb dS 



(81) 



1 's-rS> 2 



This becomes much simpler when we transform the stochastic differential equation 
to the variable x = In S using Ito's lemma to give 

(83) 



a 



dx(t) = [r- — )dt + adW(t) 



with the Hamiltonian 



1 2 ^_ / _a*\ d_ 
2° dx 2 + \ 2 j dx 



and much simpler Lagrangian 



1 /. a 



2\ 2 



(84) 



(85) 



We will consider a somewhat more general stochastic process for the Black-Scholes 
stochastic differential equation given by 

(86) 



dx(t) = ( V(x) -—)dt + adW(t) 
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which represents a situation where r is given by V(x). We can recover the Black- 
Scholes result by setting V(x) = r and a general potential V(x) is interesting as a 
mathematical exercise. 

To accomodate the discounting of all assets by the money market account, we have 
to include a killing term in the stochastic differential equation so that all expectations 
are discounted by exp(— f dtV(x(t))). This changes the Hamiltonian to 

r 2 & + H - t) I - (87) 

The reason the Hamiltonian above differs by an overall sign from the Black-Scholes 
Hamiltonian is because in pricing of options, one is considering the backward Fokker- 
Planck equation that results from the stochastic Langevin equation. 

The Hamiltonian above can be simulated by numerically integrating the Black- 
Scholes stochastic differential equation with r replaced by V(x) and putting in the 
discount factor explicitly when calculating expectations. 

Alternatively, general barrier options can be considered by keeping the Black- 
Scholes process for x but introducing an extra killing term V(x). In that case, we 
get the Hamiltonian 

*=-TS? + <T- r >J; +r+v M (88) 

where the killing term r is included because even the plain vanilla option price must 
be discounted by the money market account to get its current price. This can also be 
simulated by numerically integrating the Black-Scholes stochastic differential equa- 
tion and putting in the discount factor explicitly when calculating expectations. 
In general, therefore, option pricing with the Hamiltonian 

d =-y|? + <T- r >J; +r+vw (89) 

is equivalent to solving the Ito stochastic differential equation 

dx = (r - dt + °dW(t) (90) 

where W(t) is a standard Wiener process and then taking all expectations discounted 
by exp(— J(r + V(x(t)))dt) jHj. If we are using the Hamiltonian 

6 = -T& + iT- v Wm + v v (91) 

this becomes equivalent to solving the stochastic differential equation 

2 

dx = (V(x) - °—)dt + adW(t) (92) 
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with discounting now being done by the factor exp(— f V(x(t))dt). 

This can be done in a straightforward manner numerically. The stochastic 
differential equation can be solved by the Euler method to give sample paths 
x(ti),0 < i < N where to = 0, U = C = ie and the discounting of the final value of 
the option can be done with the factor given above. More explicitly, the simulation 
is done by numerically integrating the Ito stochastic differential equation (JHOj) using 



x i+ i = Xi + (r - ^-)(e) + a^/eZ 



(93) 



or (|9^jl using 



x 



i+ i = Xi + V(xi) 



a" 
~2 



(94) 



where Xi is x(tj) and Z are normally distributed random numbers with unit variance. 
Such random numbers can be obtained from the usual uniformly distributed random 
numbers in the range [0, 1) by the transformation 



Ci 

(2 



\J —1 ln£i cos 27r£ 2 
\/ — 2 In £i sin 27r^2 



(95) 
(96) 



where £i and £2 are distributed uniformly on [0, 1) and £i and (2 are now normally 
distributed with unit variance. 

The modified call option price desribed by ([88)1 is then given by 



E 



exp 



-tT 




X N 



-K)_ 



(97) 



where the x are generated using (}93j) . The modified call option prices described by 
(JHIIJ) are given by 



E 



N-l 



exp 



i=0 




■l-N 



-K)_ 



(98) 



where the x are generated using (|94|). 



9 Results of Numerical Simulations 

The initial stock price So is assumed to be 100. We make use of the variable x = In S 
as explained above and we set xo = In So- 

There are two ways we use the potential as explained in the previous section. In 
the first case, the killing term (coefficient of constant) in the Hamiltonian is r + V(x) 
while the drift term, the coefficient of J^, is still r — ^-. This corresponds to a constant 



20 



interest rate r and where the option payoff at expiry is defined in a path dependent 

way as 



g[S) = exp 



V{S{t))dtj S{T) - K 



(99) 



where (y) + stands for max(0, y). In this case, the potential V and r are separated as 
they have different interpretations. 

In the second case, the interest rate is assumed to be given as a function of the 
underlying security price and is equal to V(x). In this case, the Hamiltonian has V 
in both the drift and killing (coefficient of constant) terms. The option payoff g (S) 
is then path-independent and is the usual call option payoff max(0, S — K). 

Prices with potential = 1 for x<xO and for x>xO (S0=100) 

t=l year, r = 0.05, volatility = 0.25/year, 50,000 configurations, 128 time steps 

—i | i | , | , | i | , | i | i— 




100 110 120 

Strike price 



150 



Figure 4: Plot of call option price against strike price for a purely discounting poten- 
tial of form (jlUUj) and for the Black-Scholes case. 

Figures IU El and [3 refer to the potential as used in the first sense. The interest 
rate was fixed at 5%, the volatility a 2 at 0.25/year and the time to expiry at one 
year. In figure El we show a plot of the call option price versus the strike price with 



V{x) 



In figure the potential used is 



V(x) 



1 X <= Xq 
X > Xq 



(100) 




(101) 
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Prices with potential = 1 for x>xO and for x<xO (S0=100) 

t=l year, r = 0.05, volatility = 0.25/year, 50,000 configurations, 128 time steps 
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Figure 5: Plot of call option price against strike price for a purely discounting poten- 
tial of form (jlUlj) and for the Black-Scholes case. 



while in figure 13 it is 



and in figure it is 



*"(•'■) = { 1 X< X ° (102) 
U x > x 



I x <= xo 
V(x) = \ (103) 

1—1 X > Xq 

Figures |H1 El EH an d El refer to the potential as used in the second sense. The 
volatility at 0.25/year and the time to expiry at one year. In figure |H1 we show a plot 
of the call option price versus the strike price with 

V(^) = | 1 - 05 X<=X ° (104) 

[0.05 x > x 

In figure El the potential used is 

fo.05 x <= x 
V{x) = < (105) 

1.05 x > Xq 



while in figure El it is 



. -0.95 x<=x 
V[x) - ^ (106) 
0.05 x > x 
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Prices with potential = -1 for x<xO and for x>xO (S0=100) 

t=l year, r = 0.05, volatility = 0.25/year, 50,000 configurations, 128 time steps 
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Figure 6: Plot of call option price against strike price for a purely discounting poten- 
tial of form (jl(J2j) and for the Black-Scholes case. 



and in figure ^2 it is 

0.05 X <= Xq 

—0.95 x > xo 

The last two cases are unrealistic in practice as interest rates can never go negative. 



V{x) 



(107) 



10 Conclusions 

The path integral formulation of financial instruments, as shown in this work, is a 
promising approach to the pricing of derivative products which shows a remarkable 
flexibility. We have presented several applications of the method and have provided 
a general strategy -based on the use of a potential in the modeling of barrier options- 
to analyze these instruments. 

Specifically, we have compared Langevin simulations and Monte Carlo simulations 
and shown that a discount on the price of these options has indeed a rather simple 
interpretation in terms of paths of the underlying security. We believe that these 
strategies will turn out to be very effective for the simulation of complex portfolios, 
as well as for the inclusion of constraints in the evolution of these derivatives. 
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Prices with potential = -1 for x>xO and for x<xO (S0=100) 

t=l year, r = 0.05, volatility = 0.25/year, 50,000 configurations, 128 time steps 



80 



60 



40 



20 



70 



80 



90 



Price with potential 
Black-Scholes price 



100 110 120 

Strike price 



130 



140 



150 



Figure 7: Plot of call option price against strike price for a purely discounting poten- 
tial of form (jl(J3|) and for the Black-Scholes case. 
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11 Appendix A 

In this appendix we discuss some technical issues regarding the spectrum of the 
eigenvalue equation for the pricing of asian options. 

The general structure of the stationary Schrodinger equation for an Asian Option 
is of the form 



H — E) ip(x 



a 



2 d 



,<J , d 

+ ( r) — 

2 dx 2 1 2 J dx 



which we rewrite in the form 

d 



dx 2 



8 2 
- 2a- (iKe x 



dx 



[iKe x — E)ip = 



E) V = 



:i08) 



(109) 



with a = (a 2 /2 — r)/a 2 , as defined above. The velocity dependent term is elimi- 
nated as in the usual Black-Scholes model by factorizing an overall e ax term in the 
eigenfunctions. 
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Result for Potential 1.05 for x<x0 and 0.05 for x>x0 (S0=100) 



t=l year, volatility = 0.25/year, 50,000 configurations, 128 time steps 
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Figure 8: Plot of call option price against strike price for a potential of form (jl()4j) 
and for the Black-Scholes case. 



The eigenfunctions for the Asian Option can be expressed in terms of modified 
Bessel functions of the first kind I v {z) and I_ u (z), of complex argument z. Specifi- 
cally, we introduce the index function 



, . Va 2 - 2E , , 

u(E) = ^ (110) 

a 



then the solution is of the form 



Mx) = C^-^^T (1 - 2v{E)) I- 2u(E )(z) + C 2 e ax+ *^T (1 + 2v{E)) I 2u(E) ( 

^2K 



z = ze 



\z\ = 2- e x/2 (111) 



a 



;ii2) 



with C\ and C 2 arbitrary integration constant. Notice that the eigenfunctions have 
a branch cut from — oo to 0. 
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Result for Potential -0.95 for x>x0 and 0.05 for x<x0 (S0=100) 
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Figure 10: Plot of call option price against strike price for a potential of form ()106J) 
and for the Black-Scholes case. 



Result for Potential -0.95 for x<x0 and 0.05 for x>x0 (S0=100) 

t=l year, volatility = 0.25/year, 50,000 configurations, 128 time steps 
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Figure 11: Plot of call option price against strike price for a potential of form (|107|) 
and for the Black-Scholes case. 
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